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ABSTRACT 

The  role  of  halted  “baroclinic  modes”  in  the  central  equatorial  Pacific  is  analyzed.  It  is  found  that 
dominant  anomaly  signals  corresponding  to  “baroclinic  modes”  occur  in  the  upper  layer  of  the  equatorial 
Pacific,  in  a  two- and- a- half  layer  oceanic  model,  in  assimilated  results  of  a  simple  OGCM  and  in  the  ADCP 
observation  of  TAO.  A  second  “baroclinic  mode”  is  halted  in  the  central  equatorial  Pacific  corresponding 
to  a  positive  SST  anomaly  while  the  first  “baroclinic  mode”  propagates  eastwards  in  the  eastern  equatorial 
Pacific.  The  role  of  the  halted  second  “baroclinic  mode”  in  the  central  equatorial  Pacific  is  explained 
by  a  staged  ocean-atmosphere  interaction  mechanism  in  the  formation  of  El  Nino:  the  westerly  bursts 
in  boreal  winter  over  the  western  equatorial  Pacific  generate  the  halted  second  “baroclinic  mode”  in  the 
central  equatorial  Pacific,  leading  to  the  increase  of  heat  content  and  temperature  in  the  upper  layer  of  the 
central  Pacific  which  induces  the  shift  of  convection  from  over  the  western  equatorial  Pacific  to  the  central 
equatorial  Pacific;  another  wider,  westerly  anomaly  burst  is  induced  over  the  western  region  of  convection 
above  the  central  equatorial  Pacific  and  the  westerly  anomaly  burst  generates  the  first  “baroclinic  mode” 
propagating  to  the  eastern  equatorial  Pacific,  resulting  in  a  warm  event  in  the  eastern  equatorial  Pacific. 
The  mechanism  presented  in  this  paper  reveals  that  the  central  equatorial  Pacific  is  a  key  region  in  detecting 
the  possibility  of  ENSO  and,  by  analyzing  TAO  observation  data  of  ocean  currents  and  temperature  in 
the  central  equatorial  Pacific,  in  predicting  the  coming  of  an  El  Nino  several  months  ahead. 
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1.  Introduction 

The  gravest  baroclinic  waves  have  been  shown  to 
play  important  roles  in  the  ocean  adjustment  (Philan¬ 
der,  1990).  With  the  lateral  boundary  effect  consid¬ 
ered  and  coupled  to  atmosphere,  the  gravest  oceanic 
baroclinic  waves  have  been  successfully  applied  in  ex¬ 
plaining  the  dynamics  of  the  ENSO  (El  Nino/Southern 
Oscillation)  event.  ENSO  phenomena  have  been  stud¬ 
ied  worldwide  for  over  a  decade,  and  much  progress  has 
been  made  in  realizing  its  mechanism  and  sometimes 
in  giving  a  successful  prediction.  Under  the  frame  of 
large-scale  atmosphere  and  ocean  interaction,  theories 
and  hypotheses  have  been  presented  to  explain  the  de¬ 
velopment  and  quasi-cycle  of  the  ENSO  event  (e.g., 
Hirst,  1986;  Battisti,  1988).  Typical  ENSO  events 
summarized  by  Rassmuson  and  Carpenter  (1982)  were 


well  simulated  in  some  air-sea  coupled  models  (e.g.,  Ze- 
biak  and  Cane,  1987).  The  features  of  this  kind  of  El 
Nino  occurred  in  the  coupled  model,  such  as  the  west¬ 
ward  displacement  of  the  positive  SST  anomaly  and 
the  anomalous  westerly  response  to  the  SST  anomaly 
migration.  Yet  there  are  still  some  problems  that  need 
to  be  solved.  For  example,  since  1982-1983,  El  Nino 
events  always  appear  to  have  the  characteristic  of  the 
positive  SST  anomaly  moving  eastward  from  the  west 
of  the  Pacific  basin  to  the  east.  Using  an  ocean  GCM 
forced  by  the  observed  wind  field,  Philander  (1983) 
better  simulated  the  1982-1983  El  Nino  event.  Their 
simulations  gave  the  impression  that  the  El  Nino  was 
a  forced  event  but  no  exact  physical  processes  were 
given.  Figures  1  and  2  are  observed  results  from  1993 
to  1995  and  from  1996  to  1998  copied  from  the  “Cli¬ 
mate  Diagnostics  Bulletin” . 
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Fig.  1.  (a)  Zonal  wind  anomaly,  (b)  SST  anomaly  and  (c)  20°  C  isotherm  depth  anomaly  from  May  1993  to  April 
1995  (taken  from  Climate  Diagnostics  Bulletin). 
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July  1998  (taken  from  NOA A/National  Weather  Service). 


It  can  be  seen  from  these  figures  that,  before  the 
El  Nino  event  develops,  positive  SST  anomalies  first 
occur  in  the  central  equatorial  Pacific  and  then  spread 
eastwards.  In  the  atmosphere,  before  the  central  Pa¬ 
cific  warms,  westerly  anomalies  spread  over  the  equa¬ 
torial  western  Pacific.  The  timescale  for  the  westerly 
anomaly  is  intraseasonal.  Lau  and  Chan  (1986)  pre¬ 
sented  a  hypothesis  that  the  intraseasonal  timescale 
oscillations  in  the  tropical  atmosphere  might  be  a  trig¬ 


ger  for  the  El  Nino  event.  In  studying  the  1996/1997 
ENSO  event,  McPhaden  et  al.  (1999)  and  van  Old- 
enborgh  (2000)  found  that  the  onset  of  the  El  Nino 
was  linked  to  eastward  propagating  equatorial  Kelvin 
waves  forced  by  intraseasonal  atmospheric  oscillations. 
Since  the  intraseasonal  timescale  westerly  burst  hap¬ 
pens  many  times  over  the  equatorial  Pacific  each  year, 
what  mechanisms  are  there  in  the  connections  between 
the  atmospheric  intraseasonal  westerly  wind  burst  and 
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Fig.  3.  Longitude-depth  cross  section  of  long-time  mean 
water  temperature  in  April  along  the  equator  made  by 
Kessler  (taken  from  http:/ /www. pmel.noaa.gov/tao/chm/ 
clim.html). 
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Fig.  4.  Longitude-time  cross  section  of  westerly  burst 
perturbed  zonal  currents  of  two  baroclinic  modes  in  a  two- 
and-a-half  layer  tropical  ocean  model  when  the  upper  layer 
is  initially  set  at  75  m  and  the  second  layer  set  at  175  m 
(upper  layer,  positive  values  are  shaded,  units:  cms-1). 

the  occurrence  of  the  El  Nino  event?  Sun  at  al  (1997), 
Sun  and  Li  (1998)  has  analyzed  the  1994  El  Nino  event 
and  found  that  the  stationary  phenomena  of  the  pos¬ 
itive  SST  anomaly  and  sea  level  anomaly  in  the  cen¬ 
tral  equatorial  Pacific  may  be  closely  related  with  the 
El  Nino  event.  A  hypothesis  of  a  staged  atmosphere- 
ocean  interaction  mechanism  was  presented  in  their 
study.  In  this  paper,  by  analyzing  the  oceanic  re¬ 
sponse  results  from  a  two-and-a-half  layer  tropical 
ocean  model  and  assimilated  results  from  an  OGCM, 


we  found  that  the  halting  of  the  second  baroclinic 
mode  in  the  central  equatorial  Pacific  really  caused 
the  increase  of  the  SST  anomaly  in  the  boreal  winter 
in  the  central  equatorial  Pacific.  The  next  westerly 
may  be  connected  with  the  seasonal  cycle  because, 
in  the  boreal  spring,  when  the  ITCZ  is  at  its  posi¬ 
tion  near  the  equator,  convection  will  be  more  favored 
over  the  warmer  water  of  the  central  equatorial  Pacific, 
thus  leading  to  another  intraseasonal  westerly  burst 
over  the  western- central  equatorial  Pacific.  The  sec¬ 
ond  westerly  burst  will  generate  more  powerful  baro¬ 
clinic  mode  propagation  to  the  eastern  equatorial  Pa¬ 
cific,  leading  to  an  El  Nino  event.  Thus,  in  the  tropical 
atmosphere-ocean  system,  staged  interaction  processes 
are  likely  to  exist. 

2.  The  halted  mode  in  a  two-and-a-half  layer 
model  and  the  effect  of  background  wind 

Dewitte  et  al.  (1999)  studied  the  linear  vertical 
modes  near  the  equator  using  an  OGCM  and  found 
that  the  characteristics  of  the  vertical  modes  vary  more 
in  space  than  in  time.  The  first  two  modes  contribute 
with  comparable  amplitudes  but  with  different  spatial 
distributions  in  the  equatorial  wave- guide.  Boulanger 
(2001)  investigated  how  well  an  oceanic  model  whose 
dynamics  are  based  on  long  equatorial  waves  can  sim¬ 
ulate  the  large-scale  surface  zonal  current  variability 
and  validated  an  ocean  linear  model,  which  can  be  run 
with  two  or  three  layers,  against  several  sets  of  obser¬ 
vations  in  the  Pacific  Ocean.  Figure  3  is  a  longitude- 
depth  section  of  long-term  mean  water  temperature 


Longitude  (E) 

Fig.  5.  Longitude-time  cross  section  of  westerly  burst 
perturbed  zonal  currents  of  two  baroclinic  modes  in  a  two- 
and-a-half  layer  tropical  ocean  model  when  the  upper  layer 
is  initially  set  at  75  m  and  the  second  layer  set  at  175  m 
(lower  layer,  positive  values  are  shaded,  units:  cms-1). 
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Fig.  6.  Longitude-time  cross  section  of  westerly  burst 
perturbed  zonal  currents  mode  in  a  one- and- a- half  layer 
tropical  ocean  model  when  the  upper  layer  is  initially  set 
at  150  m  (positive  values  are  shaded,  units:  cms-1). 


Fig.  7.  Longitude-time  cross  section  of  westerly  burst 
perturbed  zonal  currents  of  two  baroclinic  modes  in  a 
two- and- a- half  layer  tropical  ocean  model  starting  from  a 
nearly  equilibrium  state  of  1080  days  of  integration  (upper 
layer,  positive  values  are  shaded,  units:  cms-1). 

at  May  of  2002  along  the  Eqator.  If  we  consider  the 
mixed  layer  as  the  upper  layer,  the  thermocline  layer 
as  the  second  layer  and  the  layers  below  the  thermo¬ 
cline  as  the  lower  layer,  in  the  simplest  point  of  view, 
the  tropical  ocean  can  be  divided  into  three  layers  in 
the  western  tropical  basin  and  two  layers  in  the  eastern 
tropical  basin. 

The  oceanic  part  of  the  model  in  the  atmosphere- 
ocean  system  presented  by  Zebiak  and  Cane  (1987)  is 
actually  two  layers.  Only  one  baroclinic  mode  existed 
in  their  model.  The  model  used  in  this  paper  is  a  two- 
and-a-half  layer  model  initially  created  by  McCreary 
and  Yu  (1992).  It  can  be  easily  reduced  to  a  one-and- 
a-half  layer  model  by  a  simple  change  of  parameter. 
Some  modifications  of  coefficients  and  damping  terms 
were  made  in  this  study  in  getting  the  results  of  the 
current  response  to  a  25-day  westerly  burst  in  Figs.  4 


and  5.  Wind  forcing  was  over  an  area  of  10°  longi¬ 
tude  (from  10°E  to  190°E)  and  6°  latitude  symmet¬ 
rically  across  the  equator  initially  starting  with  calm 
conditions.  It  can  be  seen  from  the  figures  that  the  in- 
traseasonal  westerly  burst  can  generate  two  baroclinic 
modes  in  this  layered  model.  Both  modes  propagate 
eastwards  but  the  second  one  propagates  much  slower 
than  the  first  mode.  The  differences  between  Figs.  4 
and  5  are  apprent.  The  current  induced  by  the  second 
baroclinic  mode  has  an  opposite  direction  in  the  upper 
layer  compared  to  the  second  layer.  When  the  model 
is  reduced  to  one  and  a  half  layers  (Fig.  6),  there  is 
only  one  baroclinic  mode.  Note  that  when  the  depth 
of  the  active  layer  is  taken  as  the  typical  depth  of  the 
thermocline  in  the  eastern  part  as  in  Fig.  3,  the  propa¬ 
gation  speed  in  the  model  ocean  is  similar  to  the  speed 
of  the  20°  C  isotherm  depth  anomaly  signal  in  the  east¬ 
ern  equatorial  Pacific  in  Figs.  1  and  2.  In  the  experi¬ 
ment  to  extract  these  baroclinic  waves  from  the  back¬ 
ground,  the  model  ocean  was  first  integrated  for  1080 
days  to  reach  the  equilibrium  state  forced  by  a  merid¬ 
ional  uniform  westward  wind  field  (McCreary  and  Yu, 
1992)  with  maximum  stress  at  61°N  (the  grid  in  the 
experiment  is  1°  x  1°).  Then,  an  anomalous  patch  of 
westerly  wind  was  added  for  25  days  to  the  region  from 
grids  10  to  19  in  a  zonal  direction  symmetrically  about 
the  equator  for  6  grids  wide  (only  the  wind  anomaly 
is  used  in  creating  Figs.  4  and  5).  It  is  pronounced 
in  Fig.  7  that  the  second  baroclinic  mode  was  halted 
in  the  central  region  of  the  model  basin.  Correspond¬ 
ing  to  the  ceasing  of  the  propagation  of  the  wave,  the 
amplitude  of  the  halted  mode  was  increased  and  the 


Longitude  (E) 

Fig.  8.  Longitude-time  cross  section  of  westerly  burst  per¬ 
turbed  temperature  anomaly  of  two  baroclinic  modes  in  a 
two-and-a-half  layer  tropical  ocean  model  starting  from  a 
nearly  equilibrium  state  of  1080  days  of  integration  (upper 
layer,  positive  values  are  shaded,  units:  °C). 
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Fig.  9.  Longitude-time  cross  section  of  westerly  burst 
perturbed  zonal  currents  of  two  baroclinic  modes  in  a 
two- and- a- half  layer  tropical  ocean  model  starting  from 
a  nearly  equilibrium  state  of  1080  days  of  integration  [(the 
background  wind  stress  is  0.85  times  stronger  than  in  Fig. 
7)  upper  layer,  positive  values  are  shaded,  unit:  cms-1]. 

temperature  of  the  upper  layer  also  increased  by  nearly 
1°C  at  its  maximum  (Fig.  8). 

When  the  strength  of  the  background  wind  was  re¬ 
duced  to  85%  of  that  used  in  generating  Fig.  7,  the 
location  of  the  baroclinic  mode  halting  was  moved  to 
the  region  near  the  eastern  boundary  (Fig.  9).  This 
reveals  that,  in  the  ocean  adjustment,  except  for  the 
first  Kelvin  baroclinic  mode,  the  other  higher  vertical 
Kelvin  mode  can  be  much  affected  by  the  background 
of  the  climate.  By  examining  the  model  experiments, 
we  found  that  the  halted  second  baroclinic  mode  in 
the  model  results  also  has  the  characteristics  of  a  pos¬ 
itive  upper  layer  depth  anomaly  corresponding  to  the 
positive  temperature  anomaly. 

In  section  3  of  this  paper,  we  examine  the  assimi¬ 
lated  results  of  vertically  averaged  SODA  data  to  dis¬ 
cuss  the  baroclinic  modes. 

3.  Analysis  of  assimilating  data  along  the  equa¬ 
torial  Pacific 

Monthly  data  assimilated  by  a  simple  ocean  data 
assimilation  (SODA)  analysis  (Carton  et  ah,  2000) 
were  used  to  analyze  the  baroclinic  modes.  The  grid 
interval  in  the  dataset  is  1°  in  the  zonal  direction  and 
variable  in  the  meridional  direction. 

The  anomalous  current  induced  by  the  second 
baroclinic  mode  has  two  zero  points  in  the  vertical 
direction.  In  the  upper  layer,  the  current  related  with 
the  second  vertical  mode  is  eastward  while  in  the  sec¬ 
ond  layer  (regarded  as  the  vertical  average  of  the  ther- 
mocline)  the  direction  is  westward.  In  the  vertical  av¬ 
erage,  the  time  serials  of  anomalous  currents,  anoma¬ 
lous  wind  stress  and  water  temperature  near  the  sur¬ 
face  were  calculated  (Figs.  10-12). 


Figure  10  gives  the  vertically- averaged  zonal  cur¬ 
rent  anomaly  at  165°E  and  the  near-surface  temper¬ 
ature  anomaly  averaged  from  170°E  to  170°W.  The 
thick  solid  line  is  for  the  current  anomaly  in  the  upper 
100  m,  the  thin  solid  line  for  the  current  anomaly  in  the 
layer  between  depths  of  100  m  to  250  m  and  the  dot¬ 
ted  line  for  the  average  water  temperature  anomaly.  It 
shows  that,  except  in  the  years  1951,  1986/1987,  1992, 
in  the  previous  boreal  winter  of  an  El  Nino  year,  there 
was  a  current  anomaly  structure  like  the  structure  of 
the  second  baroclinic  mode.  Following  with  this  type 
of  signal,  the  positive  SST  anomaly  in  the  vicinity  of 
the  data  line  is  also  increased. 

It  can  also  be  seen  from  Fig.  10  that  the  second 
baroclinic  signal  is  seasonally  dependent  because  of  the 
same  characteristics  in  wind  stress.  Liu  and  Qin  (1994) 
found  that  the  energy  of  the  intraseasonal  timescale 
oscillation  was  stronger  in  the  boreal  winter  over  the 
western  tropical  Pacific  before  the  El  Nino  event.  From 
the  study  in  this  section,  it  is  easy  to  realize  how 
strongly  the  intraseasonal  wind  stress  variation  affects 
the  El  Nino  event  by  examining  the  zonal  current  and 
SST  anomalies  in  the  central  basin  of  the  equatorial 
Pacific. 

In  Fig.  10,  there  are  15  years  with  such  character¬ 
istics.  These  years  do  not  include  the  El  Nino  years 
like  1951,  1986-87  or  1992,  but  include  1990,  a  year 
which  seemed  to  be  an  El  Nino  year  but  turned  out 
not  to  be.  In  almost  all  of  these  14  El  Nino  years, 
the  central  equatorial  Pacific  warming  anomaly  with 
the  second  baroclinic  mode  were  several  months  (1-6 
months)  ahead  of  El  Nino  event. 

Figure  11  gives  the  averaged  current  anomalies  at 
140°W  in  the  upper  layer  above  50  m,  in  the  layers 
between  50  m  and  200  m,  and  in  the  layer  below  200 
m.  Unlike  the  structure  at  165°E,  there  is  no  struc¬ 
ture  such  as  the  second  baroclinic  mode  in  the  mixed 
layer  and  thermocline.  Variations  of  current  anomalies 
in  the  upper  layer  and  the  second  layer  are  approxi¬ 
mately  in  phase  while  the  current  anomaly  in  the  lower 
layer  is  nearly  always  out  of  phase.  Thus,  at  the  lo¬ 
cation  of  0°,  140° W,  the  signal  appears  to  have  the 
structure  of  the  first  baroclinic  mode. 

Figure  12  is  similar  to  Fig.  10  except  the  wind  stress 
averaged  from  145°E  to  165°E  is  replaced  by  the  SST 
anomaly.  The  phase  of  the  wind  stress  variation  is  a 
little  ahead  of  the  phase  of  the  second  baroclinic  mode 
like  signals  during  most  El  Nino  years.  The  wind  stress 
must  be  responsible  for  the  occurrence  of  such  signals. 
It  is  implied  that  the  intraseasonal  timescale  westerly 
burst  can  generate  the  baroclinic  mode  anomaly  in  the 
central  equatorial  Pacific.  A  numerical  study  made  by 
van  Oldenborgh  (2000)  also  obtained  such  a  result  in 
the  1996/1997  El  Nino  event. 
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Fig.  10.  Vertically- averaged  zonal  current  anomaly  at 
165°E  and  near-surface  temperature  anomaly  averaged 
from  170°E  to  170°W.  The  thick  line  is  for  the  layer  from 
the  surface  to  a  depth  of  100  m,  the  thin  line  for  the  layer 
between  depths  100  m  and  250  m  and  the  dotted  line  for 
the  temperature  anomaly  (values  x20). 
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Fig.  12.  Vertically-averaged  zonal  current  anomaly  at 
165°  E  and  sea  surface  wind  stress  anomaly  averaged  over 
2°S  to  2°N,  145°E  to  165°E.  The  thick  line  is  for  the  layer 
from  the  surface  to  a  depth  of  100  m,  the  thin  line  for  the 
layer  between  depths  100  m  and  250  m  and  the  dotted  line 
for  the  wind  stress  anomaly  (values  xl50). 
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Fig.  11.  Vertically- averaged  zonal  current  anomaly  at 
140°W.  The  thick  line  is  for  the  layer  from  the  surface  to 
a  depth  of  50  m,  the  thin  line  for  the  layer  between  depths 
50  m  to  200  m  and  the  dotted  line  for  the  layer  between 
depths  200  m  and  400  m. 

The  central  equatorial  Pacific  may  be  a  key  re¬ 
gion  in  the  development  of  El  Nino  events.  Boulanger 
(2001)  found  that  the  addition  of  a  second  baroclinic 
mode  always  improves  the  simulation  of  both  the  sea 
level  and  the  surface  currents  in  an  OGCM,  especially 
in  the  central  western  Pacific. 


Fig.  13.  Vertically-averaged  zonal  current  anomaly  at 
165°  E  and  sea  surface  wind  stress  anomaly  averaged  over 
2°S  to  2°N,  160°E  to  180°E.  The  thick  line  is  for  the  layer 
from  the  surface  to  a  depth  of  100  m,  the  thin  line  for  the 
layer  between  depths  100  m  and  250  m  and  the  dotted  line 
for  the  wind  stress  anomaly  (values  xl50). 

4.  The  effect  of  the  second  baroclinic  mode 
signals  on  the  wind  stress  anomaly  over  the 
western-central  equatorial  Pacific 

When  the  regional  averaged  wind  stress  is  taken 


NO.  1 


SUN  ET  AL. 


51 


Fig.  14.  A  sketch  of  the  staged  atmosphere-ocean  inter¬ 
action  related  with  the  abnormal  warming  in  the  central 
equatorial  Pacific.  Above  and  below  the  dashed  line  are 
the  first  and  second  stages,  respectively. 


Fig.  15.  Zonal  current  anomalies  at  0°,  165° E  observed 
by  TAO.  The  thick  line  is  the  result  averaged  in  the  upper 
100  m,  and  the  thin  line  is  the  result  averaged  from  100  m 
to  250  m  (1993-1999). 

within  2°S  to  2°N,  145°E  to  165°E,  instead  of  2°S 
to  2°N,  160°E  to  180°E,  the  phase  of  the  wind  stress 
anomaly  is  changed  (Fig.  13).  In  most  El  Nino  years, 
the  anomalous  signal  of  the  oceanic  baroclinic  mode 
appears  earlier  than  the  increase  of  the  wind  stress 
during  the  boreal  winter  or  spring.  This  is  impor¬ 
tant.  It  had  long  been  pointed  out  that  the  location  of 
the  ITCZ  was  close  to  the  equator  in  the  time  leading 


up  to  an  El  Nino  event.  From  Figs.  1,  2  and  8,  the 
SST  anomaly  caused  by  the  second  baroclinic  mode 
is  around  1°C,  lasting  2  to  3  months,  corresponding 
to  the  actual  SST  beyond  28° C;  the  increased  heating 
by  warm  water  and  the  ITCZ  jointly  together  may  fa¬ 
vor  the  shifting  of  the  deep  convective  center  from  over 
the  western  tropical  Pacific  to  the  region  near  the  date 
line,  leading  to  the  prevailing  of  the  westerly  anomaly 
over  the  west  of  the  convective  region.  The  mecha¬ 
nism  discussed  here  is  the  atmosphere-ocean  interac¬ 
tion,  but  it  contains  two  different  stages:  the  tropical 
ocean  responds  to  the  intraseasonal  westerly  first  by 
a  halted  second  baroclinic  mode  and  then,  with  the 
help  of  the  ITCZ,  the  tropical  atmosphere  responds  to 
the  increased  heating  in  the  central  Pacific  by  a  second 
westerly  anomaly,  the  induced  westerly  in  turn  drives 
the  ocean  again,  generating  the  first  mode’s  propaga¬ 
tion  to  the  eastern  basin  of  the  equatorial  Pacific. 

Chen  et  al.  (2001)  studied  the  relation  between 
the  intraseasonal  oscillation  and  ENSO,  finding  that 
the  intraseasonal  oscillation  is  closely  related  with  the 
SST  in  the  Nino-3  region.  Sun  et  al.  (1997)  presented 
such  a  mechanism  (Fig.  14),  but  the  second  stage  was 
not  proved  and  was  examined  without  considering  the 
role  of  the  seasonal  position  of  the  ITCZ.  The  study 
in  this  paper  completes  the  whole  picture  and  proves 
that  the  scenario  is  possible. 

5.  Second-baroclinic-mode  signals  in  the  TAO 
observation  and  as  a  simple  prediction  index 
of  El  Nino 

Figure  15  shows  the  zonal  current  anomaly  ob¬ 
served  by  ADCP  at  165°E  on  the  equator.  It  is  also 
obvious  that  at  the  beginning  of  the  El  Nino  event, 
there  were  vertical  structures  like  the  second  baroclinic 
mode  being  observed  in  the  years  1993,  1994  and  1996. 
Thus  the  assimilated  results  are  reliable. 

In  Fig.  6,  there  are  more  instances  of  the  second 
baroclinic  mode  than  in  El  Nino  years.  For  instance, 
in  the  years  of  1999,  2000  and  2001,  there  were  sec¬ 
ond  baroclinic  modes  in  the  currents  but  the  spatially- 
averaged  surface  temperature  anomaly  around  the 
date  line  was  also  below  zero.  This  implies  that  pro¬ 
cesses  with  larger  spatial  and  temporal  scales  may  be 
other  factors  that  need  to  be  considered.  Because  the 
temperature  in  the  equatorial  region  is  more  closely  re¬ 
lated  with  the  meridional  cell  in  the  longer  time  scales 
than  the  second  baroclinic  mode,  the  interrelation  be¬ 
tween  the  halted  mode  and  the  meridional  cell  will  be 
discussed  in  a  later  paper. 

Since  there  is  an  ADCP  current  and  temperature 
observation  system  in  the  tropical  Pacific,  it  will  be 
convenient  to  judge  a  coming  El  Nino  by  examining 
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Fig.  16.  The  first  EOF  mode  of  the  orthogonal  eigenvector  (upper  panel)  and  principal  component  (lower  panel) 
at  longitudes  140° W  (left)  and  165° E  (right). 


the  ADCP  messages  from  buoy  located  at  165°E  and 
near  the  date  line. 

6.  Empirical  orthogonal  function  (EOF)  analy¬ 
sis 

The  EOF  analysis  results  for  the  zonal  current 
anomaly  and  their  energy  spectrum  at  longitudes 
165°E  and  140°W  are  presented  in  Fig.  16.  It  can 
be  seen  that  the  structure  of  the  first  dominant  mode 
is  precisely  compatible  with  the  climate  features  re¬ 
vealed  in  Fig.  3.  At  165°E,  the  thickness  of  the  positive 
zonal  current  anomaly  was  around  100-150  m,  while 


at  140°W,  the  upper  layer  thickness  became  nearly 
300  m.  As  for  their  amplitudes  of  fluctuation,  there 
was  more  vibration  energy  at  165°  E  than  at  140°W.  A 
correlation  analysis  between  150°E  and  175°E  showed 
(figure  not  presented)  that  their  simultaneous  corre¬ 
lation  was  0.58,  and  the  correlation  with  175°E  at  1- 
month  lag  was  0.47,  with  a  statistical  confidence  above 
95%.  The  correlation  between  175°E  and  150°E  at  1- 
month  lag  was  only  0.21,  a  relation  that  is  quite  poor. 
From  the  EOF  results,  it  may  be  concluded  that  the 
dynamics  of  baroclinic  waves  in  the  western  equato¬ 
rial  Pacific  are  different  from  those  in  the  eastern  Pa¬ 
cific.  More  fluctuation  energy  of  the  second  baroclinic 
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waves  is  found  in  the  western  tropical  Pacific  than  in 
the  eastern  tropical  Pacific. 

Using  Geosat  data,  Perigaud  (1990)  studied  the  sea 
level  oscillations  along  the  two  shear  fronts  of  the  Pa¬ 
cific  North  Equatorial  Countercurrent.  Filtered  sig¬ 
nals  of  50-90  days  with  the  strongest  amplitude  were 
found  east  of  160°W  in  their  study.  In  our  EOF  anal¬ 
ysis,  the  strongest  amplitudes  of  the  dominant  mode 
were  found  near  the  date  line,  where  the  significantly 
highest  frequency  of  the  variation  had  a  period  of  more 
than  3  months.  The  baroclinic  mode  analysis  in  this 
paper  seems  more  gravel  than  the  work  of  Perigaud 
(1990). 

7.  Conclusions 

The  second  baroclinic  mode  observed  in  the  cen¬ 
tral  equatorial  Pacific  might  play  an  important  role  in 
the  development  of  the  El  Nino  event.  The  intrasea- 
sonal  atmospheric  westerly  burst  can  influence  the  El 
Nino  event  through  the  staged  atmosphere-ocean  in¬ 
teraction  mechanism.  The  role  of  the  halted  second 
“baroclinic  mode”  in  the  central  equatorial  Pacific  as 
a  key  factor  in  the  formation  of  El  Nino  is  explained 
in  the  following  staged  ocean- atmosphere  interaction 
mechanism. 

The  westerly  bursts  in  boreal  winter  over  the 
western  equatorial  Pacific  generate  the  halted  second 
“baroclinic  mode”  in  the  central  equatorial  Pacific,  the 
increase  of  heat  content  and  temperature  in  the  up¬ 
per  layer  of  the  central  Pacific  induces  the  shift  of 
convection  from  over  the  western  equatorial  Pacific 
to  the  central  equatorial  Pacific,  which  leads  to  an¬ 
other  westerly  anomaly  burst  spreading  over  the  west¬ 
ern  and  central  equatorial  Pacific.  The  subsequent 
westerly  anomaly  burst  generates  the  first  “baroclinic 
mode”  propagation  to  the  eastern  equatorial  Pacific, 
which,  at  last,  results  in  an  El  Nino  event.  The  re¬ 
sults  presented  here  reveal  that  the  central  equato¬ 
rial  Pacific  is  a  key  region  in  detecting  the  possibil¬ 
ity  of  the  ENSO  and,  by  analyzing  currents  and  tem¬ 
perature  data  of  TOGA-TAO  (Tropical  Ocean  and 
Global  Atmosphere-Tropical  Atmosphere  Ocean)  ob¬ 
servations  in  the  central  equatorial  Pacific,  a  coming  El 
Nino  event  may  probably  be  predicted  several  months 
earlier. 
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